The Dynamical Cluster Approximation (DCA) versus the Cellular Dynamical Mean 
Field Theory (CDMFT) in strongly correlated electrons systems 
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We are commenting on the article Phys. Rev. B 65, 155112 (2002) by G. Biroli and G. Kotliar in 
which they make a comparison between two cluster techniques, the Cellular Dynamical Mean Field 
Theory (CDMFT) and the Dynamical Cluster Approximation (DCA). Based upon an incorrect 
implementation of the DCA technique in their work, they conclude that the CDMFT is a faster 
converging technique than the DCA. We present the correct DCA prescription for the particular 
model Hamiltonian studied in their article and conclude that the DCA, once implemented correctly, 
is a faster converging technique for the quantities averaged over the cluster. We also refer to their 
latest response to our comment where they argue that instead of averaging over the cluster, local 
observables should be calculated in the bulk of the cluster which indeed makes them converge much 
faster in the CDMFT than in the DCA. We however show that in their original work, the authors 
themselves use the cluster averaged quantities to draw their conclusions in favor of using the CDMFT 
over the DCA. 
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In their article, [lj G. Biroli and G. Kotliar compare two 
cluster methods, the Dynamical Cluster Approximation 
(DCA) and the Cellular Dynamical Mean Field The- 
ory (CDMFT) for the simplified one- dimensional large-N 
model^ and conclude that the CDMFT converges faster 
with cluster size than the DCA. 

This is a surprising result and in contradiction with exact 
results we have previously published in Ref. Q where we 
have shown that for large linear cluster sizes L c , the DCA 
converges with corrections of 0(1/1%) while the CDMFT 
converges with corrections of (D(l/L c ). The objective of 
this comment is to resolve this controversy. 
The difference in scaling behaviors of the DCA and the 
CDMFT is a direct consequence of their different bound- 
ary conditions of the cluster. The DCA 0, || and the 
CDMFT are both approximative methods incorporat- 
ing nonlocal correlations in correlated electron systems 
relinquished by fully local techniques such as the Dy- 
namical Mean Field Approximation (DMFA) HSHM 
and the Coherent Potential Approximation (CPA). [ll| 
The DCA is formulated in momentum space by dividing 
the first Brillouin zone into a number of subcells. It re- 
duces the complexity of the problem by coarse- graining 
over all the momenta within a subcell. As a result, the 
real lattice is mapped onto a cluster with periodic bound- 
ary conditions coupled to a mean field. The CDMFT on 
the other hand is formulated in real space and maps the 
lattice onto a cluster with open boundary conditions. 
We realized that the implementation of the DCA algo- 
rithm in Ref. []] is incorrect as they do not coarse-grain 
the nonlocal interaction and their analysis lacks an ac- 
curate scaling comparison in terms of the cluster size for 
these two methods. In addition, instead of presenting 
physically relevant lattice quantities they plot the corre- 
sponding cluster quantities which in the DCA in general 
are not identical to their lattice counterparts. We repeat 
the calculations of Ref. pj by implementing the DCA 
correctly and conclude that for any cluster size L c , the 
DCA results are closer to the exact solution of the stud- 
ied toy model than the CDMFT results, consistent with 



our previous findings. Q 

The exact mean field treatment of the simplified one- 
dimensional large- N model Hamiltonian (Eq. 4 in Ref. [l|) 
yields 
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f{(3E k )cosk, 



(1) 

as self-consistent equations for the dispersion Ek, with 
the number of sites L, the chemical potential fj,, the in- 
verse temperature f3 and the Fermi function f(j3Ek)- In 
the dispersion Ek, one can derive the self-energy as 



S(fe) = -2 x cosfc, 



(2) 



where the cos k factor comes from the interaction between 
nearest neighbors. 

In the DCA, the first Brillouin zone is divided into L c 
subcells with a linear size of Afc = 2ir/L c . The ir- 
reducible quantities are constructed from the coarse- 
grained propagators. Both the Green function, and 
the interaction (when it is non-local) must be coarse 
grained [fiHT^. Thus, the self-energy part in Eq.|3changes 
to its DCA counterpart 



^dca{K c 



-2 X ° CA cos(K c ) , 



where 



cos(K c ) = -y^2 cos ( K c + k) , 



(3) 



(4) 



comes from coarse-graining the nearest neighbor in- 
teraction over all the momenta k within each subcell 
around the cluster momenta K c representing the subcells. 
G. Biroli and G. Kotliar in Ref. [l| (Eq. 6 for the disper- 
sion) do not coarse-grain cos k in the nonlocal interaction 
and instead, they just substitute it by cosine. This im- 
plementation of the DCA is incorrect and overestimates 
the effects of the interaction. The correct self-consistent 
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FIG. 1: x as a function of j3 (inverse temperature) at t — 1 and 
rilling n = 0.37 (fj, w 1) for both the DCA and the CDMFT at 
L c — 2, 4, 8, 16. The exact result is represented by the solid 
circles. 

DCA equations for the dispersion are 

C K c 

^K c +fc « - 2 * C0S (^ C + *) - 2Xd Cj4 cos(^ c ) + /i. (5) 

After self-consistency is achieved, one can use the cluster 
quantity Xci CA t° compute its lattice counterpart 

xPa? t A = jJ2f^ E K c + ~k)^sk. (6) 

k 

G. Biroli and G. Kotliar plot the cluster quantity Xd, CA 
instead of the physically relevant lattice quantity XuSt A 
(Fig. 1 in Ref. [Jj). In the DCA, non-local cluster quan- 
tities like Xd, CA are n °t identical to their corresponding 
lattice quantities and thus have no physical meaning. In 
the CDMFT, since there exists one to one correspon- 
dence between the lattice quantity xf a tt MFT an d the 
cluster quantity Xci DMFT i they are identical. 
By numerically solving the self-consistent Eqs.2]and[S] 
and the one for the CDMFT (Eq. 8 in Ref. Q), we can 
make a comparison between the exact solution (cf . Eq. Q] 
in this article) and the predictions by the DCA and the 
CDMFT. 

In Fig. ^ we plot x versus j3 for t = 1 and fixed filling of 
n = 0.37 (p w 1) using the correct implementation of the 
DCA (compared to the results of Fig. 1 in Ref. 0) and 
CDMFT for L c = 2, 4, 8 and 16. It is more physical to do 
the calculations at fixed filling n instead of fixed chemi- 
cal potential /z as it was done in Ref. [lj. As illustrated, 
it is clear that the DCA results converge faster than the 
CDMFT ones to the exact solution even for small L c . 



FIG. 2: x as a function of 1/L C at /3 = 5, t = 1 and filling 
n = 0.37 (/i w 1). The insets show the convergence of the 
DCA (filled circles) and CDMFT (open squares) results to the 
exact solution (solid line) beyond L c — 16 for the CDMFT 
and L c = 32 for the DCA. The dashed lines in the insets 
represent the linear and quadratic fits in 1/L C to the CDMFT 
and DCA results respectively. 




FIG. 3: x as a function of j3 (inverse temperature) at t — 
1 and fi = 1 for the CDMFT at L c = 2,4. The second 
(filled circles) and forth (filled squares) curves are the digitized 
results extracted from Fig. 1 in Ref. Q plotted versus the 
CDMFT results obtained in this work (first and third curves) . 
The last curve (open circles) refers to \ at L c = 4 in the bulk 
(center of the cluster). 

The CDMFT results converge consistently whereas for 
small L c (L c = 8 in particular) the DCA results fluc- 
tuate around the exact solution. Nevertheless, at large 
enough L c , the DCA behaves consistently and converges 
from below. We have avoided odd cluster sizes in these 
calculations because for the DCA, the periodic boundary 
conditions on the cluster suppress Q — ir fluctuations by 
frustrating the cluster. 

In Fig. |21 we plot x as a function of the inverse cluster 
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size 1 /L c computed by the DCA and the CDMFT at fixed 
inverse temperature f3 = 5 for t = 1 and n = 0.37 (fx w 1) 
and compare the results with the exact solution (solid 
line). The DCA begins to converge quadratically and 
the CDMFT linearly in 1/L C to the exact solution beyond 
L c = 32 and L c — 10 respectively (the fits in the insets). 
The scaling behavior of both the DCA and CDMFT at 
large enough L c are consistent with our previous results. 

It is noteworthy that for this Hamiltonian in particular, 
the values of L c beyond which the scaling behaviors of the 
DCA and the CDMFT are observed are relatively high. 
E.g., in the one dimensional Falicov-Kimball model at 
half filling, these scaling behaviors are manifest beyond 

The results in both Figs. and [21 indicate the superiority 
of the DCA over the CDMFT in terms of faster conver- 
gence with the cluster size L c . The DCA also significantly 
reduces the computational CPU time compared to the 
CDMFT. In the CDMFT, because of the open boundary 
conditions on the cluster, the Green function in both real 
and momentum spaces is not diagonal. Thus, expensive 
diagonalizations of large matrices are required for large 
cluster sizes. 

In their latest response to this comment, 0] G. Biroli 
and G. Kotliar argue that in the CDMFT, whenever the 
correlation length is finite, local observables such as \ in 
this problem are best estimated in the bulk than near 
the boundaries of the cluster. In fact, their values in the 
bulk converge exponentially in the CDMFT while in the 
DCA they still converge as 1/L 2 C . Because the CDMFT 
breaks the translational invariance inside the cluster, lo- 
cal observables calculated in the bulk of the cluster differ 
from those on the boundary. Provided that the system 
is far from a transition, the sites in the center of the 
CDMFT cluster couple to the mean-field only through 
propagators which fall exponentially with distance. Lo- 
cal observables, when extracted from the center of the 
cluster instead of from their flat average, hence converge 
much faster, i.e. exponentially as opposed to l/L c . As a 
matter of fact, this point was made previously by some 
of us Q; however, for most models, it is difficult to for- 
mulate a causal self energy from quantities measured on 
the center of the cluster. Nevertheless, in their original 



work in Fig.l of Ref. [ij, G. Biroli and G. Kotliar them- 
selves perform a flat averaging of the local observable x 
and by comparing its results with their incorrect DCA 
results, they conclude that the CDMFT converges better. 
To prove that the results in Fig.l in Ref. |lj are indeed 
the flat averaged \ an( i n °t its value in the bulk, we 
plot in Fig. |3 both the flat averaged and bulk \ versus 
(3 for L c = 2,4. We have also extracted the CDMFT 
curves in Fig.l of Ref. pj for L c = 2, 4 (filled circles and 
squares respectively) using Engauge digitizing package. 
Their perfect match with the flat averaged curves and 
their deviation from the bulk result for L c = 4 best indi- 
cate that x m Fig- 1 01 R- ei - Q was fiat averaged and not 
calculated in the bulk. 

Another significant difference between these two tech- 
niques appears in the calculation of the self-energy. The 
DCA approximates the lattice self-energy by a constant 
within a DCA subcell in momentum space and as a result, 
the self-energy becomes a step function in fc-space. Thus, 
one may use an approach such as the Akima spline which 
is the smoothest curve interpolating through these step 
functions to generate an analytic form for the self-energy. 
Using the Akima spline is consistent with the assumption 
that the self-energy is a smoothly varying function in 
momentum space. In the CDMFT, since the lattice self- 
energy is approximated by the cluster self-energy in real 
space, its Fourier transform to momentum space yields 
a smooth self-energy curve. Due to the open boundary 
conditions of the cluster, the translational invariance is 
violated and the self-energy in momentum space is not 
diagonal. Nevertheless, by using the Random Phase Ap- 
proximation (RPA), |15| the translational invariance can 
be restored and the self-energy becomes a diagonal func- 
tion in momentum space. 
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